Symmetries and dynamics in an AC-driven self-assembled quantum dot lens 
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Theoretical results for a single electron in multi-level system given by a lens-shape self-assembled 
quantum dot in the presence of an intense harmonic electric field are presented. A non-perturbative 
Floquet approach is used to study the dynamical localization of the particle when going beyond 
the two-level approach by introducing the full spectral level structure. It is discussed the role of 
\^ • the different quasi-energy sidebands as the parameters of the system change. It is found that the 

contribution of different drive harmonics is controlled by fine tuning of field intensity. It is also 
shown that avoided crossings in the quasi-energy spectrum are correlated with the spectral force of 
the sidebands and dynamical state localization. 
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^ I. INTRODUCTION 

There is intense activity on the experimental and theoretical understanding of the dynamical evolution of quantum 
systems exposed to strong time-dependent external fields 0, 0] • The topic has acquired further relevance in connection 
with the practical operation of devices subjected to oscillating electrical and magnetic fields at the nanoscale. Examples 
include the shift of resonances in heterostructures as ac-fields are applied (the ac-Stark effect) Q, the behavior of 
electronic bands in spatially periodic systems 0, S B 0] > an d t ne production of currents in an ac-driven quantum 
dot Q . One important effect in these systems is the strong dynamical suppression of tunneling at suitable values of 
applied ac-ficld. The coherent destruction of tunneling that appears, known as dynamical localization in the literature, 
has been well studied in two- level systems as coming from the destructive interference introduced by the drive 0, 0] , 
whenever there is a crossing of quasi-energy levels in the spectrum. Dynamical localization has been proposed as a 
tool to control the spatial location of a particle in a two- well potential [2j , and to selectively control the tunneling in 
a multiple- well system [Io| . 

■ Typical growth conditions of semiconducting quantum dots result in dots with lens geometry, and an analysis of 
£C) ' this spatial symmetry on the electronic structure is of interest [ill [T^ |. In this work we explore the problem of 

, periodic driving force and dynamical localization in a realistic level structure that describes self-assembled quantum 
dots in semiconductors. This requires that we extend the Floquet formalism to self-assembled quantum dots with 
lens shape. We also establish the importance of incorporating the multi-level structure of a real system and identify 

■ in this complex level structure the conditions for dynamical localization. 

We find the realistic lens shape to be crucial in the description of the dynamics, as the spatial non-separability of 
the state plays a relevant role. Consideration of the multilevel structure present in typical quantum dots is a vital 
requirement for the correct description of the dynamical response of carriers. It is essential to go well beyond the 
consideration of only two active levels to fully describe the Floquet quasi-energy and the time evolution of electrons 
in realistic quantum dots. Even for weak driving forces or frequencies, the description of the dynamics requires 
the inclusion of many different states in order to achieve a fully converged description of the time evolution. The 
quantum lens geometry makes for a complex and rich theoretical description of the problem fl3| . Interestingly, we 
show that the real dot system allows the generation of higher harmonics of the driving frequency, with intensity 
that is fully dependent on the amplitude of the drive, which could be used for its generation. The phenomenon of 
dynamical localization is shown to remain for suitable values of driving field, with strongly diminished localization at 
high intensity fields. 
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II. FORMALISM AND RESULT. 



We consider a typical self-assembled quantum dot (SAQD) with lens symmetry of circular cross section of radius 
a and maximum height b which is harmonically driven by an electric field along the axial symmetry z of the lens 
with intensity F and frequency u, F = F sin Lot z. Assuming that the electron is described by an isotropic band with 
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effective mass m* , the dynamics of the system is governed by the time-dependent Schrodinger equation 

L^(r,t) = ( V 2 - eFz -rsmut -ifi—) ^(r,t) = 0, (1) 

\ 2m* at J 

where, in spherical coordinates, we have z ■ r — rcosO. The space of functions where the operatior L of Eq. ([I]) 
is defined, corresponds to those spatio-temporal functions which are bounded functions defined in the real space IZ3 
of the lens domain and are also periodic functions on time with period r = 2tt/uj. The solution of Eq. ([T]) can be 
obtained following the standard Floquet theory Q where ^(r^t) is written as 

i>(r,t)=e- t£t / h <p(r,t), (2) 

which allows to rewrite Eq. {T]) as an eigenvalue problem with the same operator L for a real-valued eigenenergy e 
(called quasi-energy) and eigenfunction (p(r,t) which fulfills the periodic condition ip(r,t) = ip(r,t + t). The periodic 
time part of the function is expanded as 

00 Qinujt 

l n \/2'K u> 
n— — 00 V / 

Furthermore, the function ifn = Q mUjt if is also a solution with quasi-energy e n = e + nhu> [J, 0] ■ These additional 
solutions, a consequence of the time periodicity, have been called replicas or sidebands. By subtracting a suitable 
integer multiple of Tilo, the quasi-energy e is mapped onto the first Brillouin zone (FBZ) which, in units of Tiw, will 
be scaled to [0,1] for convenience. As the ac-field is selected along the z-axis, the z-component of angular momentum 
with quantum number m is preserved. Then, the spatial component u n (r) in each subspace with a given value of m 

is expanded as a linear combination of the complete set of functions j for the tridimensional lens domain, as 

reported in Ref. [l2| . Finally it is obtained an infinite eigenvalue problem for the quasi-energy e and weight coefficients 
C n ,N,m- The complexity of the problem depends on the values of the field intensity F and frequency w. The solution 
is sought in terms of a truncated basis set [N = 1, N max ) which is made as large as needed in order to reach the 
desired convergence. Along this paper it is used a lens domain with ratio b/a = 0.71, 25 energy levels in the expansion 
of u n (r), and for each energy level, the index n for the replicas in Eq. (J3]) was taken as n = 0, ±1, ±2, ±20. 

Starting from an initial state at t = 0, the carrier particle is induced to explore the complete spectrum of the system 
when the ac-field is connected. How this occurs as a function of time can be analyzed by following the evolution of 
an initial electron state f (r) = f(r, t = 0), which is written as a linear combination of Floquet states 

f(r,t)= ApV P (r,t), (4) 

{P}£FBZ 

with coefficients Ap fixed by the initial conditions. The summation is taken on the First Brillouin Zone (FBZ), 
including all the replicas [l[ . At F = the label P indicates the index of quasi-energy e at,„ with nth replica and level 
N. For simplicity we consider that the initial state has well-defined z-component of the angular momentum in such 
a way that only states with the same m are considered. Then, Eq. can be cast in the following way 



/ m (r,t)=^A Ar (t)<: ) (r), (5) 

N 



where 



00 p i(nujt - 7r/2) . ,„ w , t 

A N (t)= £ f; A P C n , N , m (P) e e -te(P)t/h (6) 



{P}eFBZ n=-oo 



^/2tt/uj 



The function Ajv(i) contains all the dynamical information due to the presence of the ac electric field. Thus, P/v(i) = 

|Ajv(i)| 2 gives, at a given time t, the probability of finding the system in the state „ • The typical definition of 
dynamical localization considers that the carrier does not evolve away from its initial spatial state or configuration, 
as the effective tunneling amplitude from site to site is suppressed [2|. Here, we monitor dynamical localization by 
the corresponding quantum probability of finding the particle in its initial state as a function of time. In what follows 

we assume that the initial state f (f) is the zero-field ground state ®^m( r ) with N = 1 and m = 0. 
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FIG. 1: Oscillation Amplitude as a function of the dimensionless field intensity eFa/Tiu) for b/a = 0.71 and Q, = 100. Insets 
show the probability Pi of finding the particle in the initial state as a function of time for different values of eFa/hui. 
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FIG. 2: Same as in Fig. [Qfor = 300. 

In order to quantify the degree of localization, we study the Oscillation Amplitude (OA) as a function of the field 
intensity, where the OA is defined as OA = P±(t = 0) —Pmini where P m in is the minimum value that P\(t) takes 
over a long interval of time (100 time units in our case) [14j. In Fig. [1] and [5] the OA is plotted as a function of the 
dimensionless parameter eFa/Tku where Q — Tilu/E is a dimensionless frequency with E Q — h 2 /(2m*a 2 ) as a unit of 
energy. At zero field the OA is naturally always zero and increases with the field strength, since the latter induces 
mixing of the states and it forces the system to explore larger regions of the eigenvalue spectrum of the system away 
from the initial state. The OA reaches the maximum value (OA = 1) at certain value of the reduced field intensity 
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FIG. 3: a) Small window of the FBZ for the quasi-energies as a function of the reduced field intensity, b) Intensity sidebands 
pp(n) for P — 2 and P = 13 with n — 0. c) Idem as b) but now with n — —3. Calculations are done for b/a = 0.71 and 

n = ioo. 



and for eFa/hcu > 11 it is near unity, indicating that the probability of the system remaining in the initial state is 
zero. However, the OA strongly decreases at some particular field intensities (for example at eFa/hcu = 4 in Fig. [1] 
and at eFa/hcu = 7.2 and 9.0 in Fig. [2]), meaning that Pi(t) never goes to zero at those values of the field intensity, 
and quasi-localization of the system can be identified. At the same time and to study explicitly the behavior of Pi 
on time, different panels inside Fig. Q] and [2] show the time evolution P\{t) during one hundred time units r = H/Eq 
at different values of the dimensionless field intensity. According to Eq. (JSJ), the wave function is a multi-periodic 
function and presents strong oscillating behavior as seen in the inset of the figures. Notice, moreover, that at small 
value of reduced frequency (f2 = 100), the mixture of the spectrum is strong for eFa/hcu = 3.6 in Fig. Q] and for 
eFa/hcu — 4 a condition for quasi-localization is reached. Instead, for higher values of the frequency (f2 = 300), a 
lower slope for the increase of the OA is obtained as seen in Fig. [2J It is also reported two values for quasi-localization 
at eFa/hcu — 7.2 and 9.0. Notice also that the stronger the intensity the lower the values of the probability, a 
consequence of the strong level mixing as discussed previously. 

On the other hand, the interaction among different lens levels in the quasi-energies as function of the field amplitude 
can be analyzed in more detail using the sideband intensity or spectral force pp, which amounts to the weight of the 
n-th sidebands on the tpp(r,t) Floquet state. Thus, according to the expansion of u n (r), the spectral force is given 
by 

pp(n) = ]T|auv, m (P)| 2 . (7) 

N 

Notice that at F = 0, pp(n) = pnin). On the other hand, for F > 0, the spectral force of a given P indicates 
the weight of the different lens states in the full Floquet expansion . Figure [3] a) shows a small window of the 
quasi-energy spectrum, where some anticrossings are present at O = 100. The system parameters are the same as 
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in Fig. Q] In Figs. [3]b) and[3]c) the intensity sidebands pp(n) for the quasi-energies P — 2 and P = 13 have been 
plotted, with n = in Figs. [3]b) and n = —3 in Figs. [3]c) respectively. We observe in the figures the effects in the 
spectral force p p (n) caused by the the anticrossing between quasi-energies P = 2 and P = 13 around eFa/hcu ~ 1.8. 
For values up to eFa/fko ~ 1.8, pp=2{n — 0) and pp=i3(n = —3) show a large intensity as expected due to its 
corresponding F — limit, while pp = \^{n = 0) and pp=2{n = —3) are barely noticeable near F ~ 0. However, at the 
corresponding anticrossing, a strong mixing of the states takes place, and pp = \z{n = 0) and pp—2(n = —3) increase 
rapidly while pp=2{n = 0) and pp=i3(n — —3) strongly decrease for eFa/fko > 1.8, i.e., their corresponding strengths 
are inverted after the anticrossing. Thus, the spectral weights for n = and n = —3 in this case, and in general for 
all the spectral contributions, are exchanged between the Floquet states at the anticrossing. Most importantly, as 
the field increases, so does the interlevel mixture and the weights of the various replicas become nearly identical, so 
that the amplitude of their contribution are similar. We emphasize that this nearly homogeneous distribution of the 
spectral force over many different replicas results in substantially different time evolution of the driven system. We 
then anticipate that time averages of physical observable could be substantially different than when only the lowest 
two lens levels are considered 1151 . 



III. CONCLUSIONS 



We have analyzed the time evolution of an electron in self-assembled quantum dots with lens shape in the presence 
of intense radiation along to the rotational axis of the lens. We exploit this axial symmetry of the lens domain to 
solve this complex time-dependent problem. This realistic driven single-electron system has been studied over a wide 
range of field amplitudes. We have demonstrated that consideration of the typical two- level approximation yields an 
incomplete description even at moderate fields and frequencies. We have also calculated the complex quasi-energy 
spectra that result in this problem, and analyzed the anticrossings that appear in terms of the interaction among 
zero-field states and their replicas. We find that these anticrossings are associated with strong shifts in the spectral 
weights for the Floquet states between two quasi-energies, and that for larger field intensities the spectral weights are 
distributed homogeneously among a wide range of sidebands. This strong dependence indicates that the appearance 
of different drive harmonics in the response of the system could be easily controlled by the field strength. It is 
interesting to consider the possibility of utilizing such lens shape quantum dots in strong ac-fields as a source of 
different harmonics. 

In order to study whether dynamical localization prevails under these much more complicated conditions of multi- 
level dynamics, we have studied the time evolution of the system prepared initially in the zero-field ground state. 
We find that at field intensities for which a quasi-energy anticrossing appears in the spectrum, some degree of 
state localization is observed; in those circumstances, the probability of finding the system in the initial state never 
goes completely to zero, but it reaches a minimal value. This incomplete dynamical localization has been analyzed 
quantitatively by means of the oscillation amplitude for different system parameters. Similar to the case of two-level 
systems, we find that the dynamical localization, already precarious at high frequency, disappears for lower frequency 
values. The analysis and results we present are important for the interpretation of experimental data, and suggest 
further theoretical work to assess the relevance of multi-level structures in realistic systems. 
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